Scalar Fields as Dark Matter in Spiral Galaxies 
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We present a model for the dark matter in spiral galaxies, which is a result of a static and 
axial symmetric exact solution of the Einstein-Dilaton theory. We suposse that dark matter is a 
scalar field endowed with a scalar potential. We obtain that a) the effective energy density goes like 
l/(r 2 + rl) and b) the resulting circular velocity profile of tests particles is in good agreement with 
the observed one. 
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PACS numbers: 95.30.S, 04.50, 95.35 

One of the greatest puzzles of physics at the moment 
is without doubt the existence of dark matter in cosmos. 
The experimental fact that the galaxy masses measured 
with dynamical methods do not coincide with their lu- 
minous galaxy masses gives rise to the existence of a 
great amount of dark matter in galaxies, galaxy clusters 
and superclusters. At the present time, cosmological 
observations indicate that the universe is filled out with 
about 90 percent of dark matter, whose nature till now 
remains unexplained. Recently some authors have pro- 
posed the scalar field as a candidate for dark matter 
in cosmos in some sense the inflationary cosmo- 

logical model proposes the scalar field as cosmological 
dark matter as well. These models consider scalar-tensor 
theories of gravity where one is able to add mass terms 
to the total energy density of the space-time. All mod- 
ern unifying field theories also contain scalar fields. For 
example, scalar fields are fundamental fields in Kaluza- 
Klein and Superstring theories, because such fields ap- 
pear in a natural way after dimensional reduction. In 
both theories the scalar field could be endowed with an 
exponential scalar potential HI]], in particular, when 
we deal with 5-dimensional Kaluza-Klein theories, the 
Lagrangian density reads £5 = R5 + A being A a 5- 
dimensional cosmological constant. After dimensional 
reduction and a conformal transformation one obtains 
the density £4 = -Ri + 2(V0) 2 + e" 2 /^A, where <f> is 
the scalar field which actually states that an exponen- 
tial potential appears in a natural way in this theory. 
An analogous procedure stablishes that in the low en- 
ergy limit of Superstring theory one gets a similar result 
1^,1). fn general one obtains the Lagrangian from high- 
dimensional theories £ 4 = -i? 4 + 2(V$) 2 + e" 2Q<f, A, 
therefore here we will restrict our selves to an exponen- 
tial scalar potential. In this letter we show a possible 



model for the dark matter in spiral galaxies, supposing 
that such matter is of scalar nature. 

There is a common approach to explain the rotation 
curves in spiral galaxies called Modified Newtonian Dy- 
namics (MOND) |^,||], which basically consists of modify- 
ing the Newton's law of attraction for small accelerations 
by adding terms to the gravitational potential. In this 
way, by adjusting some free parameters for each galaxy, 
one can reproduce the asymptotic behavior of the rota- 
tion curves. However it appears to be artificial because 
it is nothing but a mere correction of Newton's law, we 
are unable to know neither where the parameters and 
the correction terms come from, nor why nature behaves 
like that and therefore which is the Newton's law at cos- 
mological scale for instance. 

A convincing phenomenological model for galactic dark 
matter is the called Isothermal Halo Model (IHM), which 
assumes the dark matter to be a self-gravitating ball of 
ideal gas (made of any kind of particles) at a uniform 
temperature kT = \mdrnVc-, - being rrid m the mass of 
each particle and v c its velocity - which eventually pro- 
duces a dark matter distribution going as ~ I/r 2 , im- 
plying in this way an increasing mass M(r) ~ r. Then, 
by assuming that a galaxy is a system in equilibrium 
(GM/r 2 — v 2 /r) the velocity of particles surrounding 
the profile above should produce flat rotation curves into 
a region where the dark matter dominates, i.e. at large 
radii when one considers as exponential distribution of 
luminous matter as usual Qj. 

Observational data show that galaxies are composed 
by almost 90% of dark matter (6[j^]. This is so because 
the kinematics inside the dark matter dominated region 
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is not consistent with the predictions of Newtonian the- 
ory, which explains well the dynamics of the luminous 
sector of the galaxy but predicts a keplerian falling off 
for the rotation curve. The region of the galaxy we are 
interested in is that in which the dark matter determines 
the kinematics of test particles. So we can suppose that 
luminous matter does not contribute in a very important 
way to the total energy density of the matter that deter- 
mines the bahavior of particles in the mentioned region, 
instead the scalar matter will be the main contributor 
to it. Thus, as a first approximation we can neglect the 
baryonic matter contribution to the total energy density 
for the explanation of assimptotic rotation curves. 

On the other hand, the exact symmetry of the dark 
halo stills unknown, but it is reasonable to suppose that 
it is symmetric with respect to the rotation axis of the 
galaxy. In this letter we let the symmetry of the halo as 
general as we can, so we choose it to be axial symmetric. 
Furthermore, the rotation of the galaxy does not affect 
the motion of test particles around the galaxy, dragging 
effects in the halo of the galaxy should be too small to af- 
fect the tests particles (stars) traveling around the galaxy. 
Hence, in the region of our interest we can suppose the 
space-time to be static, given that the circular velocity 
of stars (like the sun) of about 230 Km/s seems not to 
be affected by the rotation of the galaxy and we can con- 
sider a time reversal symmetry of the space-time. So, the 
model we are dealing with will be given by the gravita- 
tional interaction modified by a scalar field and a scalar 
potential. The model consists of the following action 

Sb- + 2(V$) 2 -Vm (1) 
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which could be the four-dimensional action for the 
Kaluza-Klcin or the Low Energy Superstings theory 
without electromagnetic field, and where we have added 
a term which contains the scalar potential. In this action 
R is the scalar curvature, $ is the scalar field, ko — 16 J G 
and \/—g is the determinant of the metric. The most gen- 
eral static and axial symmetric line element compatible 
with this action, written in the Papapetrou form is 

I 



ds z 



'2k 



f 



(dzdz) + W 2 d(/) 2 } - f c 2 dt 2 



(2) 



where z := p + and z := p — i£ and the functions 
/, W and k depend only on p and £. This metric rep- 
resents the symmetries posted above. The application 
of the variational principle to ([!]) gives rise to the field 
equations 



(3) 



which are the Klein-Gordon and Einstein's field equations 
respectively; /i, v = 0,1,2,3. Using the harmonic maps 



ansatz we find the following Poisson like structure 

for the above equations |12|| 
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and a similar equation for k z , with z instead of z, where 
A is the Laplace operator such that for any function 
h = h{z,z): Ah := (Wh }Z ), z + {Wh )S ) tZ . Moreover, 
A = ln(/) is interpreted as the gravitational potential. 

If one assumes that A and $ depend only on W(z,z), 
the set of equations (0) appears in a more tractable form 
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and the corresponding expression for k Zl where now the 
operator D means Dh{W) = Wh" + 2h' V7i = h(W), 
and ' denotes derivative with respect to W. Equations 
(^- ||) constitute a system of coupled differential equa- 
tions because y/—g = We 2k ~ x /2. However it is evident 
that once we have expressions for A and k can be 
integrated by quadratures. Moreover, using the third of 
these equations A and $ obey differential equations with 
W being the independent variable. 

In order to find an exact solution, we substitute 
KQy/—gV(<&) from (|t]) into ([|) and (^), (remember that 
4^ = —2aV) and obtain two decoupled differential equa- 
tions, one for A and another for $. We solve these two 
differential equations and substitute the solution into (|^). 
We thus find that a solution of the system (|[]|) is given 
by 

A = ln(M) + ln(/ ) 
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where /o and $o are integration constants and W = M 
is a function restricted only by the condition 
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whose solutions are M = Z(z)Z(z), where Z is an ar- 
bitrary function. The reader can check that (^) is a 
solution of the field equations substituting the set (^) 
into (^) using the metric (j^). 

In what follows we study the circular trajectories of 
a test particle on the equatorial plane taking the space- 
time (g) as the background. The motion equation of a 
test particle in the space-time (Q) can be derived from 
the Lagrangian 
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dr 
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This Lagrangian contains two constants of motion, the 
angular momentum per unit of mass 



W 2 d(f> 
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and the total energy per unit of mass of the test particle 



f c 2 — - A 
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where r is the proper time of the test particle. An 
observer falling freely into the galaxy, with coordinates 
(p, £, 4>, t), will have a line element given by 



ds 2 
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e 2h (p 2 + ( 2 ) + W 2 ^ 2 }- f}c 2 dt 2 
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The velocity v a — (p, is the three- velocity of the 

test particle, where a dot means derivative with respect 
to t, the time measured by the free falling observer. The 
squared velocity v 2 is then 

^ = ffabW V = ^(p 2 + C 2 ) + ^ 2 , (15) 

where a, 6 = 1,2,3. Using ( |l4| ) into ( p"3| ) we obtain an 
expression for the squared energy 



A = 
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We are interested in test particles (stars) moving on the 
equatorial plane £ = and the equation of motion de- 
rived from the geodesies of metric (p|) reads 



1 2fc / dp 
f {dr 
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where we have used the conservation equations (|l2j) and 
(|l3|) . Equation ( |l7| ) determines the trajectory of a test 
particle around the equator of the galaxy, in this tra- 
jectory A and B remain constant. If we change of test 
particle, we could have another constants of motion A 
and B determining the trajectory of the new particle. A 
spiral galaxy is practically a disc of stars traveling around 
the equatorial plane of the galaxy in circular trajectories 
in the period of observation, altough it had to be formed 
from enormous clouds of gas going around a symmetry 
axis with average values of A and B. Thus for a circular 
trajectory p = 0, the equation of motion transforms into 



B 2 f 
W 2 



c 2 f 



(18) 



This last equation determines the circular trajectories of 
test particles travelling on the equator of the galaxy. Us- 
ing ( |l8| ) and ( |l6| ) we find an expression for B in terms of 
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since v 2 <C c 2 . Now using (|T^) one concludes that for our 
solution (g) v 2 = f§B 2 , i.e. 



vdm = foB, 



(20) 



where we call v — * vdm the contribution of our dark 
matter to the circular velocity of a star. 

When Z = z our solution in Boyer-Lindquist coordi- 
nates p — v 'r 2 — 2ar + b 2 sin 9, £ = (r — a) cos 6 reads 

(l-f) 2 + ^, dr 2 



ds 2 
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where K 2 = b 2 — a 2 and r$ only scales. The effective 
energy density pdm of (||) is given by the expression 
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and plays the role of our dark matter density profile. 

Keeping in mind that this is only the contribution of 
dark matter to the energy density we are in conditions 
to compare these results with those given by measure- 
ments. In order to do so we recall the paper by Begeman 
et al. where an energy density profile of the IHM 
^(r) = p r 2 /(r 2 + r%) for dark matter is used, being r c a 
core radius. It is evident that this profile is a particular 
case of the expression we present here, namely, for mat- 
ter localized on the equator of the galaxy. So, we can 
fit some of the free parameters of metric ([H]) comparing 
these two profiles. We set b — r c , a = and ^° r " — por 2 . 




FIG. 1. The circular velocity profiles of four spiral galax- 
ies. Solid lines represent the total circular velocity (vc), 
long-dashed is the contribution of the dark matter to the total 
velocity (vdm) and the short-dashed curves is the contribu- 
tion of luminous matter (vl); finally the dots represent the 
observational data. The units are in (Km/s) in the verti- 
cal-axis and in (Kpc) in the horizontal-axis. 

Let us model the circular velocity profile due to the 
luminous matter of the disc in a spiral galaxy by the 
function 

1 97x 122 

v 2 = u 2 (iiopt ) / 3_____ (23) 

which is the aproximatted model for the Universal Ro- 
tation Curves (URC) as was propossed by Persic et al. 
H for an exponential thin disc, valid for a sample of 967 
galaxies; in this expression x = r/R opt , the parameter 
(3 — VL(R pt)/v(Ropt) being R opt the radius into which 
it is contained the 83% of the onservable mass of the 
galaxy and v the observed circular velocity. 

We can suppose that luminous matter near the center 
of a galaxy behaves like in Newtonian mechanics. Thus 
with the luminous velocity ( p3| ) it is now easy to calcu- 
late the angular momentum (per unity of mass) of the 
test particle in the luminous matter dominated region 

B = v L D, (24) 



where D is the distance between the center of the 
galaxy and the test particle. For our metric, D = 
J ds, keeping 9, <f> and t constant, we obtain D = 
y/ (r 2 — 2ar + b 2 )/ f r Q . Observe that after we have de- 
termined the dark matter energy density pdm, B is 
uniquely determined by vl via (pi]); it is easy to show 
that B in (|24]) equals that of equations ( |Tt| - |20| ) by in- 
cluding a luminous newtonian component into the radial 
geodesic equation |l3| . Therefore ( f20| ) and ( pij ) imply the 
total circular velocity 

vc = \fl + v 2 DM =v LX /l + /op - 2ar + b 2 ) (25) 

expression that should fit the observed rotation curves. 
In order to do so, we present in Fig. 1 the plots contain- 
ing the fittings of four spiral galaxies and into Table | 
the values of the parameters /o and b keeping a = and 
the scale tq — 1. From it we see that the agreement of 
the resulting circular velocity profiles given by the scalar 
field as dark matter and the observed is very good not 
only far away from the center of the galaxy but inside the 
part of the galaxy dominated by luminous matter as well. 

The criterion used to choose the sample was the ra- 
tio of the dark to luminous mass inside r25, which was 
selected to be ~ 1 in order to test our model by using 
" dark enough galaxies" . The plots shown in Fig. 1 would 
not be enough to state that our model works, it is nec- 
essary to be consistent with the phenomenological URC 
approach, i.e. the contribution of our dark matter should 
be the same as that propossed into the URC frame which 
is strongly luminosity dependent. A formula consistent 
with (^3|) is given by || 

KrcDM = V 2 (R op t)(l - /3)(1 + l 2 )^-^ (26) 

being (3 — 0.72 + 0.44 log L/L* the same parameter as in 
(||) and 7 = l^L/L,) 1 / 5 . 

According to d20| ) and ( p4| ) the contribution of our dark 
matter is 

1 Q7 T 1.22 

vIm = fo(r 2 + b 2 )v 2 (R opt )P {x2 + Q 782)1A3 (27) 

which after using the fitting parameters of Table | the 
comparisson of both approaches are compared in Fig. 
2 from which it can be concluded that our dark matter 
model respects the luminous matter model we have used. 



NGC1560 NGC2403 




4 



FIG. 2. The contribution of dark matter to the circular ve- 
locity of test particles is shown. Solid lines emerge from the 
model described in this letter, dashed ones correspond to the 
URC approach. The discrepancies are of 8.3, 11, 23.3 and 5.2 
percent respectively. 

Some remarks can be drawn. The energy density ( p^ ) 
coincides with that required for a galaxy to explain the 
rotation curves of test particles in its halo, but in our 
model, this energy density is product of the scalar field 
and the scalar field potential, that is, this dark matter is 
produced by a $ particle. So we have shown that there is 
an exact solution which describes the rotation curves of 
particles in a spiral galaxy. The crucial point for having 
the circular velocity vdm = foB is that / ~ W in the 
solution (^). But this fact remains unaltered after con- 
formal transformations in the metric ds 2 = A(<&)ds 2 , so 
that the circular velocity vdm remains the same for all 
theories and frames related with metric (||) by conformal 
transformations. 

How does our model looks like into the cosmological 
context? When a density profile for galactic dark matter 
goes as the inverse of r 2 and it is supossed that the halo 
of a galaxy ends in the region where those of neighboring 
galaxies start, the integrated amount of galactic dark 
matter is close to that needed for the Universe to be flat 
for the observed average distance between them [Q, flat- 
ness inferred from the cosmic background radiation Q] 
and thus permiting our model to be inside the bounds. 
In fact, we have developed a cosmological model that 
considers the same theory as here (|lj) with the same 
scalar potential [^5) , which has been able to explain the 
redshifts of Type la supernovea, and all the parameters 
for structure formation lie into the ranges imposscd by 
observations Jl^Jl^j , which make us to put forward the 
model presented in this letter. 
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TABLE I. Values of the fitting parameters fo and b. Here 
also are shown the values of the quantities used. 



Galaxy 


fo 


b 


Ropt 


P 




(Kpc- 1 ) 


(Kpc) 


(Kpc) 




NGC1560 


0.0726 


2.119 


4.6 


0.344 


NGC2403 


0.0171 


5.399 


6.7 


0.546 


NGC3198 


0.0038 


12.88 


11 


0.547 


NGC6503 


0.0290 


3.035 


3.8 


0.702 
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